Nonlinear conductance in molecular devices: molecular length dependence 
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We theoretically study the electronic transport in the monolayer of dithiolated phenylene vinylene 
oligomeres coupled to the (111) surfaces of gold electrodes. We use non-equilibrium Green functions 
(NEGF) and density functional theory(DFT) implemented in the TranSIESTA package to obtain a 
full ab initio self-consistent description of the transport current through the molecular nanostructure 
with different electrochemical bias potentials. The calculated current-voltage characteristics (IVC) of 
the systems for the same contact geometry have shown a systematic decrease of the conductivity with 
the increased length of the molecules. We analyze the results in terms of transmission eigenchannels 
and find that besides the delocalization of molecular orbitals the distance between gold electrodes 
also determines the transport properties. 
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I. INTRODUCTION 

Since the pioneering work on electron transport 
through molecules by Aviram and Ratner^ there has been 
remarkable progress in the field of molecular electronics. 
With the advances of experimental tools and techniques, 
various molecular devices have shown a number of inter- 
esting effects: negative differential resistance;^ switch- 
ing properties;^ memory effects^ A number of measure- 
ments have been made where the electrical current has 
been driven through monolayer molecular films between 
metallic electrodes, including break junction^'^ and evap- 
orated electrode^ experiments. The connection of a sin- 
gle molecule to conducting electrodes has been achieved 
using STMi and break junction^ techniques. Owing to 
their stability the most common systems considered were 
thiol-ended organic molecules on gold electrodes. 

The current-voltage measurement by Reed et of 
dithiol benzene molecule (dtb) onto gold electrodes in 
break-junction experiments showed an apparent gap of 
about 0.7 V. It raised the question of its nature: Whether 
it was a Coulomb gap or it reflected the mismatch be- 
tween the contact Fermi level and the lowest unoccupied 
molecular orbital. It was viewed as a question of bonding 
between the molecule and the electrodes. 

In general, in the nonbonding, and even weak-bonding 
situations one expects no coherent transport, but in- 
stead a sequential charging of the molecule and Coulomb 
blockade effects. Recently Kubatkin et al.^ have demon- 
strated a Coulomb blockade and different charge states 
in the junction with weakly coupled (physisorbed) 0PV5 
molecule onto gold contacts. 



In the strong-bonding CclSG, Oil coherent transport is 
dominant and the staircaselike structure could appear 
owing to the succession of scattering resonances entering 
the bias window.^ In many junctions no apparent gap was 
noticed. Instead, linear current-voltage curves were re- 
ported at low bias for strongly bonded thiolate molecules. 
There are a number of possible transport mechanisms 
that give rise to linear IVC. In the evaporation electrode 
technique, for example, grainy surfaces can form with 
the likely protrusion of electrode atoms in between the 
molecules of the monolayer, thus spoiling the surface or- 
der and the dominant conduction mechanism appears to 
be hopping with the characteristic energy scale of 10-100 
meV.i° 

In order to explain experimental IVCs, the dithiol- 
benzene molecule coupled to Au(lll) surfaces was a kind 
of model system for theoretical studies. After the first 
semiempirical investigations ^ii^ ab initio calculations 
within the framework of density functional theory were 
performed. The conductance was found to be more than 
an order of magnitude larger than in experiment de- 
spite the assumption of unequal bonding of the molecule 
to the electrode surfacesii^ 

In the theoretical model of molecular junction with 
parallel perfect electrode surfaces, it came natural that 
the electronic structure of the system provided a fi- 
nite density of states at the Fermi level of the cou- 
pled system, ruling out the appearance of the gap in 
the I-V curves. The theoretically obtained conduc- 
tances were as a rule more than an order of magnitude 
larger compared with the measured ones. The differ- 
ence was related to the varying strengths of the sulfur- 
gold bonds and to the adsorption sites. Interestingly, 
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for longer molecules, quantitative agreement of theoreti- 
cal considerationsiSiiLiS is considerably improved when 
compared with the measurements. ^'^^ 

The transmission properties of the systems are de- 
termined by the electronic structure of the combined 
molecule and electrode systems. The distance between 
the electrodes has an important role in the conduction of 
the junction. In simple terms, the metal-molcculc-metal 
system can be viewed as a finite quantum well that sets 
up in between the metal surfaces. The direct tunneling 
between the electrodes, despite the possible electrode- 
induced gap states in the metal-air-metal junctions, is 
considered to be too small to contribute considerably to 
the current of the typical size molecular junctions."' 

With the molecule bridging the interelectrode gap, a 
possibility for the molecule-induced interelectrode gap 
states, or resonances appears. Those states will have 
an important impact on the transport properties. The 
induced gap states, being molecule mediated, depend on 
the length of the molecule, and on the interelectrode sep- 
aration. The transport is determined by scattering states 
which are Bloch-like states in the electrodes and molec- 
ular orbitals like states in between the electrodes. 

Previous studies of the electronic transport of the junc- 
tions with the molecules of different length included a 
number of oligomeres in a search for a good molecu- 
lar wirei^2*2i For this purpose, one has to consider the 
molecules with the same basic unit that repeats succes- 
sively with the length of the molecule. The electronic 
structure at the surface and the adsorption site geom- 
etry for that reason should be similar for all the ad- 
sorbed molecules considered. Magoga and JoachimSS 
used a semi-empirical approach in their studies of the 
tunnel transport regime of certain families of oligomeres. 
They considered the transmission at zero bias voltage 
and found that the conductance follows a simple expo- 
nential decrease with the length of the molecule. How- 
ever, in order to get a full account of the details of the 
electronic structure of the electrodes, the electrode sur- 
faces, the molecule, and the molecule-induced gap states 
on the transport properties of the junction a full ab ini- 
tio approach is needed. Moreover, as the junction is sub- 
ject to the non-equilibrium conditions the transmission 
should be obtained self-consistcntly for each bias voltage 
applied across the junction. 

In this paper, we therefore use the non-equilibrium 
Green functions technique and the density functional the- 
ory to obtain a full ab initio self-consistent description 
of the transport properties. We investigate the conduc- 
tance of the sulfur-ended phenylene vynelene oligomeres 
coupled to gold electrodes; OPVn systems, with n equal 
to the number of benzene rings in the molecule. The 
number of benzene rings of these oligomeres determines 
the length of the molecule. Owing to a strong bonding of 
sulfur on metal surfaces, particularly gold, all other pa- 
rameters of the system can be kept fixed for the systems 
with different molecules of the series. The difference in 
the transport properties of the coupled systems then re- 



flects the electrode distance in addition to the electronic 
structure of the molecule. The results will show two dif- 
ferent transport regimes^^; the one where the electrodes 
are close enough that the induced gap states give an ap- 
preciable contribution to the current, and the other where 
the separation between the electrodes is large enough to 
result in their negligible contribution. 

In Sec. ^ of the paper we first report on the model 
used and the generalities of the computational procedure. 
In Sec, mil we calculate the current through the system 
as a function of applied bias voltages and analyze them 
in terms of the evolution of transmission spectra. We 
discuss the results in Sec. II VI 



II. CALCULATIONS 

A typical molecular electronic system consists of a 
molecule coupled to two electrodes with different electro- 
chemical potentials. Depending on the experimental 
method of building, the device may vary the contact ge- 
ometry with respect to the adsorption site as well as to 
the adsorption strengthiS4 We consider an idealized ge- 
ometry with a molecule symmetrically coupled to the two 
identical gold electrodes. 



A. Computational method 

To perform the first-principle quantum modeling of 
the electronic structure under non-equilibrium condi- 
tions and to calculate the current-voltage (I-V) charac- 
teristics of the system, we used non-equilibrium Green 
functions technique based on density functional the- 
ory as implemented in the recently developed TranSI- 
ESTA simulation packagei^^ Core electrons were mod- 
eled with Troullier-Martinsi2& soft norm-conserving pseu- 
dopotentials and the valence electrons were expanded in 
a SIESTA localized basis setm^ The system was divided 
into three regions: the left and right electrodes and the 
central region. The central part contained the portion of 
physical electrodes, where all screening effects took place. 
The electronic occupations of the system was determined 
by the electrochemical potentials of the electrodes. The 
charge distribution in the electrodes corresponded to the 
bulk phases of the same material. 

The density matrix of the system under the external 
bias was calculated self-consistently within density func- 
tional theory (DFT) using the only major approxima- 
tion of the method in the choice of the exchange corre- 
lation functional. We used the local density approxima- 
tion (LDA). The virtue of the method is that it gives 
the Hamiltonian in the same form as in the empirical 
tight-binding approach making its techniques straight- 
forwardly applicable. The nonlinear current through the 
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FIG. 1: A phenylene vinylene oligomere, 0PV4, connected 
to two Au(lll) surfaces via thiolate bonds. 

contact is calculated using the Landauer formula™ 

I(Vb)^Go r" T{E,Vb)dE, (1) 

where Gq = 2e^//i is the quantum unit of conductance 
and fJ-L/R SLie electrochemical potentials of the left and 
right electrodes. With the applied bias potential eVJ,, the 
left and right potentials become /xl(V6) = ij,l{0) + eVb/2 
and fJ.R{Vb) = Mi?,(0) — 614/2, respectively. The energy 
region of the transmission spectrum that contributes to 
the current integral in Eq.l. we refer to as the bias 
window. The total transmission probability T(E, Vb) = 
^^^-^ T„(i?, Vb) for electrons incident at an energy E 
through the device under the potential bias Vb is com- 
posed of all available conduction channels with the indi- 
vidual transmission T„. 



B. OPVn molecules bonded to gold electrodes 

The aim of this paper is to investigate the electrical 
and structural properties of the monolayer of phenylene 
vinylene oligomere molecules (OPVn) coupled to gold 
electrodes via the sulfur atoms on both ends. There are 
good reasons to believe that the result for the mono- 
layer of OPVn molecules applies to the single molecule, 
too.^ 

It is usually assumed that sulfur forms thiolate bonds 
with the gold surface, which is not necessarily always 
the case. Depending on the absorption technique, if the 
strong SH bond is not cleaved and if the molecule forms a 
thiol bond with gold, a preferential adsorption site might 
be at the top with the nonzero tilt anglei^ In our cal- 
culation, however, we consider thiolate bonds at a 3-fold 
hollow site and the tilt angle of the molecular axis to the 
surface normal to be zero. We have also made the calcula- 
tion with molecule adsorbed in the on-top position on the 
gold atom at both electrodes, and found that the adsorp- 
tion at the 3-fold hollow site is energetically favorable to 



the on-top absorption by 2.6 eV. In the experimental sit- 
uation, some molecules in the self-assembled monolayer 
will be atypical bonded to various kinds of defects on 
the surface. We, however, consider an idealized situation 
and assume further that all molecules align with parallel 
phenyl rings. There is a good reason to believe that this 
is energetically more favorable, since it gives a larger tt 
overlap. 

The geometry of the system was determined as follows. 
We used a (3x3) unit cell of the Au(lll) surface for the 
electrodes and assumed that the molecule of cis confor- 
mation was chemisorbed to the surface of the electrodes, 
as shown in Fig^ The 0PV3 molecule was positioned 
symmetrically in the hollow site of both Au(lll) sur- 
faces. After positioning the molecule perpendicularly in 
the z-direction so that it has a favorable Au-S bonding 
distance, we let its molecular coordinate to relax until 
the average force on the molecular atom was less than 
0.04 eV/A. The gold atoms were kept fixed at their bulk 
positions. 

The same procedure was followed for 0PV4 and 0PV5 
molecules. The obtained distances for Au-S atoms were 
2.39A for all systems considered within the computa- 
tional accuracy, the convergence in total energy tolerance 
was O.OOOleV and the force tolerance was 0.04 eV/A. In 
order to get an additional check on the stability of the 
results with respect to the Au-S distance, we made a 
calculation for the 0PV3 system with an Au-S distance 
enlarged by 0.1 A. The minor effects on the transmis- 
sion with no consequences on the current through the 
system were obtained. The lack of the force dependence 
on the transport properties of the system with chemi- 
cally bonded contacts was observed by Cui et alii^ on 
octaneditiol molecules. 

In the TranSIESTA procedure the charge on the 
molecule is not fixed. It adjusts itself to minimize the 
free energy as the electrochemical potentials of the elec- 
trodes are changed. For the systems considered, we have 
found that molecules preserve charge neutrality to within 
0.15 e with the change of the bias voltage. 

The important aspect of the calculation is the robust- 
ness of the results to computational details. From the 
calculation by Stokbro et al.^^ for the dithiolate benzene 
molecule on Au(lll) surfaces (dtb) in the hollow posi- 
tion we know that the use of the generalized gradient 
approximation (GGA) instead of the LDA gives almost 
negligible effects on the transmission spectra. The unim- 
portance of the relaxation of the first two gold layers for 
final results was also demonstrated. As the binding ge- 
ometry is mostly determined in the first two layers of 
gold atoms, the conclusion is valid for the systems with 
0PV3, 0PV4 and 0PV5 molecules as well. 
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FIG. 2: Current and differential conductance of OPVn sys- 
tems as a function of bias. A self-consistent calculation has 
been employed for each bias voltage. 



III. TRANSMISSION PROPERTIES AND I-V 
CHARACTERISTICS 

In order to obtain I-V spectra and differential conduc- 
tance, we performed self-consistent calculations for a bias 
voltage in the range from -3.0 V to 3.0 V. Since our elec- 
tronic structure in the contact region is fully symmetric, 
the corresponding I-V spectra are symmetric with respect 
to the reversed bias. 

The general trend of I-V curves is shown in Fig[21 One 
observes that all curves have a similar general shape. 
They show a slow increase of the current with the bias 
voltage up to a certain bias and then a steep increase at a 
point where transmission resonances come into alignment 
with the bias window. 

Two characteric features of the systems can be ob- 
served. The first one is that the current in the ohmic re- 
gion of the transport is lower as the length of the molecule 
is larger. The second one is that a steep increase of the 
current of the longer molecule happens at smaller bias 
than for the short molecule. 

When compared to the results for shorter DTB 
systemic the absolute value of the current of the 0PV3 
is almost an order of magnitude smaller. It is also in 
good agreement with the current of the similar type of 
molecules iSiii 



The zero bias conductance G = jfr{^L/ii,0) is also 
much smaller than for DTB. They are 2.1 fiS for 0PV3, 
1.0 fiS for 0PV4, 0.2 fiS for 0PV5, in comparison with 
35.0 ^iS for DTB. 

In the following we analyze the obtained IVCs for 
OPV3-5 in terms of the dependence of the transition 
spectra on molecular length (number of benzene rings), 
bias voltage, and interelectrode separation. 



A. OPV3 

Chemisorption of the molecule on the gold surfaces 
leads to the change of the energy levels in both systems. 
In particular, some of the molecular levels broaden into 
a continuum. The eigenstates of the combined system 
become the scattering states which are Bloch- wavelike in 
the electrode and molecular orbital-like in the moleculei^^ 
The more delocalized the state is, more spreading the 
corresponding transmission peak shows. 




-1 1 
E-Ef (eV) 

FIG. 3; Transmission amplitude of DTB and 0PV3 
molecules on Au(lll) surfaces at zero bias voltage. MPSH 
eigenvalues are marked with circles for dtb and with triangles 
for 0PV3. 

In FigEl we compare the transmission spectra of the 
coupled dtb system with the 0PV3 system at zero bias. 
In the figure we also mark the eigenvalues of the projected 
self-consistent Hamiltonian (MPSH) onto the molecular 
orbitals. The transmission spectra are proportional to 
the density of states of the combined system. The MPSH 
levels are normally broadened and shifted in energy by 
the interaction with the gold surface. Some of the states 
have a large weight on the S atom and form a strong 
bond with the gold electrode. The others may have a 
little overlap with Au states and form a narrow peak in 
the density of states and in the transmission spectra. 

In the DTB system the huge transmission resonance 
below the Fermi level is due to a superposition of strongly 
hybridized levels, as demonstrated in Ref. 15. 

The 0PV3 system does not exhibit such a huge reso- 
nance as DTB does and it has a considerably less spectral 
weight at the Fermi level. The 0PV3 resonance below 
the Fermi level is built on the HOMO molecular state 
at the energy E — — 0.799ey. The reason for being less 
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FIG. 4: Transmission amplitude and MPSH eigenvalues of 
OPV4 and OPV5 systems at zero bias. 

pronounced than in DTB is the larger separation of Au 
electrodes in 0PV3, resulting in the smaller overall over- 
lap of gold states with the molecular states. This leads to 
smaller hybridization and less broadened levels resulting 
in lower density of states around the Fermi level. 

When the bias Vb is applied to the system the chem- 
ical potential of the electrodes changes. The calculated 
transmission spectra are plotted for a set of bias values in 
FigEl The main effect of the increased bias is that more 
of the resonance gets into the bias window. The position 
of the HOMO related transmission resonance does not 
change with respect to the average electrochemical po- 
tential as the bias is increased, but its weigth diminishes. 
The resonance gives a lower contribution to the current 
tranport at higher bias, even when the maximum of the 
resonance enters the bias window. That peculiar bias 
dependence of the HOMO resonance leads to a step-like 
behavior of the current and a dip in the conductance at 
the bias of 1.25^ as seen in Fig|21 

The transmission peaks corresponding to LUMO, and 
higher levels change linearly with the applied bias, con- 
trary to the HOMO peak. They are slowly approaching 
the electrochemical potential of the left electrode leading 
to a smaller HOMO-LUMO separation. The net result 
is the lower transmission in the HOMO-LUMO region. 
The current is considerably smaller when compared to 
the current of the smaller dtb system, particularly at low 
bias. Eventually, at high enough bias, at above approxi- 
mately 2V the LUMO resonance enters significantly the 
bias window and contributes dominantly to the current. 
This is the region of the steep increase of the current as 
seen in Fig|21 Nevertheless the current in the 0PV3 sys- 
tem is almost an order of magnitude smaller than in the 
DTB system in the region we considered. 
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FIG. 5: Transmission T(E, Vt,) of the OPV3 system as a 
function of the bias voltage (shifted vertically by 0.8 for vis- 
ibility). Energies are relative to the average electrochemical 
potential (/^l + ^ir)/2. The bias windows are marked with 
dashed lines. The amplitude corresponding to the HOMO 
levels (marked with arrows) decrease with the increased bias. 



B. OPV4 

There is a clear separation of HOMO and LUMO levels 
in the 0PV4 system for zero bias, as shown in Fig0| 
The appearance of the lower transmission of the HOMO 
related resonance and the lower density of states at Ef 
is more pronounced than in 0PV3 system. The position 
of the HOMO level is 0.685ey below the Fermi level. 
The LUMO level is I.136ey above the Fermi level and 
is less broaden giving the lower tail at the Fermi level. 
The behavior is maintained at the finite bias, as seen in 
FigEl The MPSH levels are narrower in general, and 
the resulting transmission at low bias is smaller when 
compared with the 0PV3 system as seen in Fig|21 Like 
in the 0PV3 a dip in the conductance appears at 1.25U. 



C. OPV5 

The tendency of the appearance of a sharper well- 
defined peaks in the system with the longer molecule 
is clearly demonstrated in the 0PV5 system in Fig0| 
The contribution of the underlying states to the reso- 
nance is obviously changed with respect to the 0PV4. 
The overall overlap of the molecular states with the gold 
electrode states is lower leading to the narrowing of the 
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FIG. 6: Transmission T{E, 14) of the 0PV4 system for dif- 
ferent bias voltages The notation is the same as in Figl^ 
The HOMO amplitude decreases with increased bias and the 
LUMO amplitude moves towards the average electrochemical 
potential of the electrodes. 



transmission resonances. The HOMO and LUMO reso- 
nances at zero bias are at 0.684ey below and 1.063ey 
above the Fermi level, respectively. With the increased 
bias, the HOMO transmission almost disappears, while 
the LUMO slightly moves towards the average chemical 
potential of the electrodes. The dip in the conductance 
at 1.25V, which was present in 0PV3 and 0PV4, is miss- 
ing in 0PV5 owing to a very little spectral weigth of the 
HOMO resonance. 

There is a new feature in the spectra related to the 
narrowing of the levels. An almost pure HOMO-LUMO 
gap sets up at any bias. The bias dependence of the gap is 
reflected only through its slight shift with respect to the 
average chemical potential of the electrodes, as seen in 
FigtZI The low zero-bias conductance is 0.2 fiS, and the 
current is distinctly lower compared with the preceding 
member of the OPVn series. 



IV. DISCUSSION 

The current-voltage dependence of the systems can 
be obviously understood in terms of the transmission 
resonances. The 0PV3 system at zero bias does not 
show a huge transmission resonance below the Fermi 
level, as was characteristic of a smaller DTB system. 
The tendency of diminishing transmission and density 
of state around the Fermi level is more pronounced for 
the 0PV4 system and even more for the 0PV5. Loosely 
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FIG. 7: Transmission T{E, Vt) of the 0PV5 system. The 
notation is the same as in Fig|^The behavior of the HOMO 
and LUMO amplitudes is similar to that of the OPV3 and 
OPV4 systems. Note a low density in the HOMO-LUMO 
region. 



speaking, for a short molecule there will be a strong 
molecule-mediated overlap of electrode states leading to 
wide transmission resonances as opposed to the weaker 
overlap, and more narrow transmission peaks for longer 
molecules. 

To understand this consideration better, let us con- 
sider the transmission related to the interelectrode sepa- 
ration. We define it as the distance between the surface 
planes of the left and right electrodes. In our calcula- 
tions the shortest distance is in the DTB system with 
the value 1.05 nm. Our input radius of the gold atoms 
is 3 A, and the direct overlap of the gold atom states 
from the left and the right electrodes cannot take place. 
The pure gold levels, which perfectly follow the electrode 
chemical potential, can contribute directly to the tun- 
neling process through the electrode-induced resonant 
states or even quantum well states related to gold surface 
statesi^ The properties of these states are dependent on 
the separation between the electrodes. As reported by 
Kergueris et al.,^ the tunneling conductance in the gold- 
air-gold junction is linear and of the order 0.02 at 
the barrier height of IV for the interelectrode spacing 
of 1 nm. That appears to be several orders of mag- 
nitude smaller than the typical conductance of metal- 
molecule-metal junctions. However, with the molecule 
chemisorbed in between the electrodes, the electrode in- 
duced gap states and resonances play an important role 
in transport properties. They are, of course, a part of 
the hybridized molecule-electrode states, caused by the 
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TABLE I: Interelectrode separation and zero-bias conduc- 
tance of OPVn systems, d^iec. is the distance between the sur- 
face (111) gold planes of the left and right electrodes. G{^l/b., 
0) is the conductance at the Fermi level for zero bias. 





DTB 


OPV3 


OPV4 


OPV5 




1.05 nm 


2.3 nm 


2.97 nm 


3.64 nm 


G{lXL/R, 0) 


35 /iS 


2.1 fiS 


1.0 fiS 


0.2 fiS 



molecule-mediated overlap of electrode states. We can 
tentatively call them molecule-induced gap states. 

If the electrodes are close, they contribute to the finite 
density of states and to the enhanced transmission in 
the HOMO-LUMO region of the combined systems. As 
the electrode separation is larger, the molecule-induced 
gap states loose the weight, and their contribution to the 
conducting channels of the system becomes smaller. The 
transmission peaks become narrower and the transmis- 
sion in the HOMO-LUMO region consequently smaller. 

The DTB molecule is the shortest one we have consid- 
ered. The increase of the interelectrode separation in the 
OPVn series from 2.3 nm for 0PV3 to 3.64 nm for 0PV5 
is followed by the decrease of the zero-bias conductance, 
as seen in Table 2^'^ . 

Indeed, as discussed in Section Hill the corresponding 
transmission in the HOMO-LUMO region of those junc- 
tions systematically drops down with the distance be- 
tween the electrodes, as seen in Figs|31and0] The true 
discrete nature of molecular resonances starts to govern 
the transport. The level spacing of the resonances is 
illustrated in terms of the eigenvalues of the projected 
self-consistent Hamiltonian onto the molecular orbitals 
(MPSH) in Fig|Sl The levels are more closely spaced and 
the broadening of the levels goes down with the length 
of the molecule. 

The semiconducting property with a narrow band gap 
of the TT conjugated system is well known. The gap is 
due to a distinct alternation of the bonds, with every 
second bond having double-bond character. In the poly 
p-phenylene vinylene molecule, the HOMO-LUMO gap 
of 2.5eV has been reported)^ The reduction of the gap 
obtained in this paper, as seen in FiglSJ is expected for 
the thiolate-ended molecule on gold surfaces. 

Another characteristic feature of the OPVn series is its 
zero bias conductance. One finds a remarkable difference 
between the DTB and 0PV5 systems. For the DTB, it 
is 35 ^S, while for the 0PV5 of the same contact geom- 
etry it is only 0.2/iS'. Obviously, the transmission in the 
HOMO-LUMO region appreciably drops down in the se- 
ries as seen from Fig l5l7l For an even longer molecule, 
one expects a negligible density of states around Ef and 
a vanishing transmission in the HOMO-LUMO region for 
low bias voltages. 

These two features, the fixed HOMO-LUMO separa- 
tion and the gap in the density of states around Ef, show 
that a kind of bias threshold sets up. A current through 



dtb opv3 opv4 opv5 




FIG. 8: Eigenvalues of the projected self-consistent Hamilto- 
nian of the Au(lll)-molecule-Au(lll) system onto the molec- 
ular orbitals at zero-bias. Full lines connect the HOMO and 
LUMO levels, respectively. 



the junction will start at a bias above some threshold 
value, at which the wings of the HOMO and LUMO peaks 
enter the bias window. This is a kind of bias-activated 
coherent transport. 

One expects the gap in the I-V curve to be more pro- 
nounced in the systems with weaker chemisorptions than 
in the sulfur-ended molecule on gold electrodes. 

Owing to the bias gap, a long molecule, even though 
chemisorbed on the electrodes, may show a behaviour 
of a small capacitoriSiiSi^ For nonequal bonding to 
the electrodes, this quite often happens in real junc- 
tions. This opens the possibility for a kind of sequen- 
tial transport in tt conjugated systems and charging of 
the strongly chemisorbed molecule in gate-controlled ex- 
periments. Kubatkin et al^ have demonstrated gate- 
controlled Coulomb blockade involving different charge 
states in the 0PV5 molecule of tran conformation, but 
with the protected thiol end-groups of their molecule, 
thus preventing strong chemisorption. 

In the view of this, a long molecule coupled to the 
electrodes may also serve as a testing ground for the 
boundary region between mesoscopic and nanostructure 
systems. 

We also point out that the calculations we have per- 
formed showed that the hollow position is the minimum- 
energy absorption site for all molecules in the series. 
We have taken the distance of sulfur to gold atoms as 
5% larger and found that it does not change the con- 
duction noticeablyi^ Besides the known mechanical and 
temperature stability of thiolate-ended organic molecules 
bonded to gold surfaceSji^ this robustness of conduction 
is an additional desirable property in the view of appli- 
cability in molecular electronics. 



V. CONCLUDING REMARKS 

Using the NEGF and DFT we have obtained the non- 
linear current-voltage characteristics of thiolate-ended 
phenylene vinylene oligomers coupled to Au(lll) elec- 
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trodes. In the low-bias region we have found that the 
current scales down inversely to the interelectrode sep- 
aration, with the strong tendency of the appearance of 
the I-V gap for long molecules. The reason for the gap 
formation is the narrowing of the HOMO and LUMO 
related transmission peaks with the interelectrode sep- 
aration. The nonexistence of the I-V gap in the junc- 
tions of short molecules appears to be a consequence of 
the ideal Au(lll) semi-infinite surfaces and the molecule- 
mediated gap states encountered in the calculations re- 
sulting in wide resonances and finite density of states 
around the Fermi level. The measured conductance of 
short molecules is usually smaller than the calculated 
one, which is possibly due to such ideal conditions sel- 
dom, if ever, encountered in real junctions. For longer 
molecule, the effects appear to be less efficient, and as ex- 
pected the discrepancies between the measured and cal- 



culated values are less pronounced. We point out that the 
full account of such effects could only be followed with 
the full ab initio self-consistent calculations we performed 
at each bias voltage. 

We have also found that the conductance is robust to 
the sulfur gold distance for the molecules chemisorbed in 
the local energy minimum, hollow site on gold surfaces. 
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